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ABSTRACT 

Using cosmological hydrodynamic simulations of the LCDM model, we 
present a comparison between the simulation sample and real data sample of H I 
and He II Lya transmitted flux in the absorption spectra of the QSO HE2347- 
4342. The ACDM model is successful in simultaneously explaining the statistical 
features of both H I and He II Lya transmitted flux. It includes: 1.) the power 
spectra of the transmitted flux of H I and He II can be well fitted on all scales 
> 0.28 h~^ Mpc for H, and > 1.1 h~^ Mpc for He; 2.) the Doppler parameters 
of absorption features of He II and H I are found to be turbulent-broadening; 3.) 
the ratio of He II to H I optical depths are substantially scattered, due to the 
significant effect of noise. A large part of the r^-scatter is due to the noise in the 
He II flux. However, the real data contain more low-// events than simulation 
sample. This discrepancy may indicate that the mechanism leading extra fluctu- 
ations upon the simulation data, such as a fluctuating UV radiation background, 
is needed. Yet, models of these extra fluctuations should satisfy the constraints: 
1.) if the fluctuations are Gaussian, they should be limited by the power spectra 
of observed H I and He II flux; 2.) if the fluctuations are non-Gaussian, they 
should be limited by the observed non-Gaussian features of the H I and He II 
flux. 

Subject headings: cosmology: theory - large-scale structure of the universe 
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1. Introduction 

Lya forest lines and transmitted flux of QSOs' absorption spectra provide the most 
valuable samples in studying the physical state of the intergalactic medium (IGM) and 
gravitational clustering at high redshifts. High resolution samples of QSOs' Lja absorption 
spectra are important to test models of cosmic structure formation on small scales. Recently, 
the H I Lya transmitted flux of HE2347-4342 has been used to compare with hydrodynamical 
simulation samples of the ACDM model (Jamkhedkar et al. 2005, here after Paper I). The 
results suggest that the ACDM model is successful in explaining the power spectrum and 
intermittency of the HE2347-4342 sample. There is no discrepancy between the simulated 
and observed flux fields with regards to their statistical behavior from the second to the 
eighth orders and till the comoving scales as small as about 0.28 h^^ Mpc. This result seems 
not to support the necessity of reducing the power of density perturbations relative to the 
standard ACDM model on small scales, as was imphed by the lack of dense cores in the halo's 
center given by the so-called universal density profile (Flores & Primack 1994; Swaters et al. 
2003; McGaugh et al. 2003; Zentner & Bullock 2003; Simon et al. 2003) 

The IGM is also traced by He II Lya absorption. Because the ionizing threshold of He II 
is high (54.4 eV), and recombination rate of He III is also high, the He II Lya absorption 
of IGM generally is much stronger than H I Lya. Therefore, it is expected that the He II 
Lya forest and transmitted fiux of high redshift QSOs can play a similar role as H I forests 
in constraining cosmological models, and can even yield stronger constraints on the models 
than H I forests (Zhang et al. 1995; Croft et al. 1997). However, due to the lack of He II 
data, the comparison of He II spectra between model predictions and observations could not 
be made in a similar way as for H I Lya absorption spectra. Thanks to the FU SE data of 
HE2347-4342, we can obtain moderate resolution spectra of He II Lya transmitted flux in 
the redshift range 2.0 < 2; < 2.9 (Kriss et al. 2001). It provides the possibility of making a 
similar analysis as for the H I transmitted flux. 

We cannot simply repeat the analysis as Paper I, because of a new problem: the ratio 
between the optical depths of He II and H I, 77 = 4rHeii/THi. If assuming 1.) the effect of 
thermal broadening and peculiar velocities of IGM are neghgible, and 2.) He II and H I are 
in photoionization equilibrium, the ratio 77 should basically be constant, and have a low level 
of scatter because 77 is weakly dependent on the temperature of IGM. Observations reveal, 
however, that 77 is signiflcantly scattered from pixel to pixel, or from line to line (Kriss et al. 
2001; Smette et al. 2002). For the FUSE data of HE2347-4342, the scatter of r] is from 1 to 
a few hundreds (ShuU et al. 2004), and even as high as ~ 10^ (Zheng et al. 2004). 

The large 77-scatter has been used as an indicator of the inhomogeneity of UV radiation 
background caused by the radiation transfer (RT) of the UV photons in a nonuniform density 
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field. A 3-dimensional radiation transfer calculation on the shadowing, self-shielding and 

filtering predicted that the mean of 7] should be as large as (rj) > 200 (Maselli & Ferrara 
2005). However, the observed result is (r/) < 100 (Zheng et al. 2004; Shull et al. 2004). 
Therefore, the RT effect explanation of t] scatter is far from settling. This result motivated 
us to reconsider the assumptions above mentioned. If above-mentioned assumptions 1.) and 
2.) are not hold, the ry-scatter would not be a direct measurement of the inhomogeneity of 
UV background. 

In this paper, we will study the effect of thermal broadening and peculiar velocities on 
the H I and He II Lya transmitted ffux. It has been shown recently that in nonlinear regime 
the velocity field v{x) of cosmic baryon matter consists of strong shocks on various spatial 
scales and in high and low mass density area (Kim et al. 2005). The statistical behavior 
of the velocity field is similar to a fully developed turbulence (He et al. 2006). Therefore, 
the effect of thermal broadening and pecuhar velocities of IGM would not be neghgible. We 
should, at least, estimate the imprints of the non-trivial velocity field on the Lya forests of 
He II and H I. 

The paper is organized as follows. Section 2 describes the observed data of HE2347- 
4342. Section 3 presents the method to simulate the Lya forests of He II and H I. As in 

Paper I, we use the WIGEON method of cosmological hydrodynamic simulations to produce 
the simulation samples, as this code is especially effective in capturing singular and complex 
structures (Feng et al. 2004). On the other hand. He II would be formed in hotter gas, 
such as shock heated IGM (Cen & Ostriker 1999; Dave et al. 2001), a method of effectively 
simulating shocks would be important. Velocity field effect on absorption width is shown in 
Section 4. Section 5 presents the analysis of the ratio between the optical depths of He II 
and H I Lya transmissions. The power spectrum of the H I and He II flux will be discussed 
in §6. Discussions and conclusions are given in §7. 



2. Data of HE2347-4342 

The data of H I Lya transmitted fiux used in this paper is the same as in Paper 1. 
The FUSE data of the He II Lya transmitted flux of HE2347-4342 is described in Zheng 
et al. (2004). The wavelength region is 904 — 1188 A, which corresponds to the redshift 
range 2.0 < 2; < 2.9, as the wavelength 303.78 A of He II Lya in rest frame. The spectrum 
has a constant bin size of AA = 0.025 A. In terms of the local velocity, the resolution is 
dv ~ 8.3 — 6.3 km s^^, and mean rff ~ 7 km s^^. The mean S/N is 2.14. Following the 
approach of Shull et al. (2004), we bin the data into AA = 0.05 A to reduce uncertainties 
in the effective spectrograph resolution and oversampling effect. The distance between N 
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pixels in the units of the local velocity scale is given by 5v — 2c[l — exp(— A'"di>/2c)] km s ^, 
corresponding to comoving scale D — 5v{l + z)/H{z). 

The flux in 2729 pixels, i.e., 24% of the total pixels, are less than zero. Obviously, 
the points with negative flux is unphysical, it should be excluded in the statistics below. 
For sample deleting all the pixels with negative flux, the mean transmission about 0.4, or 
effective optical depth ^ 0.9. The optical depths of He II over the ranges from 0.1 to 2.3 
are with 10% uncertainties. A better statistical measurement of the fluctuations of flux is 
given by the ratio between the optical depths of He II, thcIi, and H I, thi, for each pixel. 
The distribution of r] is scattered in the range from 0.1 to about 500, while the mean of r] is 
~ 80 (ShuU et al. 2004). 

We also assume that all absorption in the FUSE spectrum is due to He II, although 
it is subject to metal-line contaminations. Generally identifled metal-lines are connected 
with a Lyman- limit system (Smette et al. 2002). The Doppler width of metal lines are 
generally narrow with 5v < 20 km s~^. In this paper, we restrict our analysis only to scales 
5v > 30 km s~^ where metal-hne contaminations is low (Hu et al. 1995; Boksenberg et al. 
2003; Kim et al. 2004). 

3. Hydrodynamic simulation sample 
3.1. Method 

We use the cosmological hydrodynamic simulation samples produced by the Weno for 
Intergalactic medium and Galaxy Evolution and formatiON (WIGEON) code developed by 
Feng et al. (2004). It is a hybrid hydrodynamic/A^-body simulation, consisting of the WENO 
algorithm (Jiang & Shu 1996) for baryonic fluid, and A^-body simulation for particles of dark 
matter. The baryon fluid obeys the Navier-Stokes equation, and is gravitationally coupled 
with coUisionless dark matter. We have assumed a standard ACDM model, which is specified 
by the matter density parameter fim = 0.27, baryonic matter density parameter Qb = 0.044, 
cosmological constant Qa — 0.73, Hubble constant h — 0.71, the mass fluctuation erg = 0.84, 
and scale-free spectrum index n — 1. The ratio of specific heats of the IGM is 7 = 5/3. The 
transfer function is calculated using CMBFAST (Seljak & Zaldarriaga 1996). 

The simulation was performed in a periodic, cubic box of size 50 h~^M.pc with a 512^ 
grid and an equal number of dark matter particles. It starts at redshift 2; = 99. A uniform 
UV-background of ionizing photons is switched on at 2; = 6 to heat the gas and reionize the 
universe. To mimic the enhancement of temperature due to radiation transfer effects (Abel 
& Haehnelt 1999), a thermal energy of gas with T — 2 x 10^ K is added in the total energy 
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at 2; = 6. The clumpy universe would reprocess the photon spectrum from ionizing sources. 

The reprocessing UV spectrum has been calculated by Haardt & Madau (1996). We use an 
ionizing background model including QSOs and galaxies with 10% ionizing photons escape 
fraction (kindly provided by F. Haardt). Ai z = 2.5, such an ionizing background produces 
the transmission flux of H I and He II similar to observation and an average 77 ~ 72, which 
is very close to observed value. 

The atomic processes in the plasma of hydrogen and helium of primordial composition, 
including ionization, radiative cooling and heating, and the fraction of H I and He II are 
calculated in the same way as Thcuns ct al. (1998). That is, under the "optically thin" 
approximation, once density and temperature of baryon gas are given, the ionizing state of 
H and He is directly determined from the ionization-equilibrium equation. 



3.2. Samples of Lya transmitted flux 

For given sample of the fields of density, temperature and velocity of the baryon matter, 
the optical depth of H I or He II can be produced by a convolution with Voigt profile as 
follows (Bi et al. 1995; Zhang ct al. 1997) 

r.(.) =a.c|dxn.(x);^=J^V[^^^ (1) 

where i =H I or He II, (jj is the absorption cross section of Lya line, ni{x) the number 
density, 6z is the rcdshift difference between z and x, v{x) the peculiar velocity in unit 
c and hj = {2kT /niiC^Y^'^ the thermal velocity. V is the Voigt profile, which is normal- 
ized / dx{l/ ^/^Tb) V[ h(^ilz) + ^] ~ 1- The Hubble constant at redshift z is H{z) = 
Ho^/^m{^ + -2)^ + ^A- E<1-(1) shows, when the terms of thermal broadening and peculiar 
velocity are not negligible, the Lja transmission fiux depends on mass density field nni, ?^Heii 
as well as the fields of temperature and velocity of IGM. 

We produced 100 mock samples of H I and He II Lya transmitted flux at z = 2.5 with 
randomly selected lines of sight. Each mock spectrum is sampled using 2^° pixels with the 
same spectral resolution as the observation. As the corresponding comoving scale for 2^° 
pixels is larger than the simulation box size, we replicate the sample periodically. We add 
Gaussian noise to H I sample with signal-to-noise ratio, S/N=50, while He II with S/N=3. 
Figure 1 shows typical samples of H I and He II Lyct transmitted flux flelds. 
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4. Line width of He II and H I Lya absorptions 

To demonstrate the importance of the velocity field, we consider the line widths of Lja 
absorption. If the Lja absorption lines are purely thermal broadened, the line width, or 
Doppler parameter, 6(He II) of He II should be less than 6(H I) by a factor of 2. However, 
Zheng et al. (2004) found 6(Hc II) = ^6(H I) and ^ ~ 1. They concluded that the velocity 
field in IGM is dominated by turbulence. 

This result can be explained with considering the velocity field in eq.(l). A fiux field 
given by eq.(l) generally is not given by a superposition of lines with Gaussian profile, 
because the peculiar velocity v{x) is a random field. Therefore, the line width given by 
Gaussian profile fitting is not always equal to the thermal broadening b^. For instance, 
when bj is small, the Voigt profile will picks up only pixels with 5z + {l + z)v{x) ~ 0. In this 
case, the line width can be estimated by (1 + 2;)((Af )^)^/^, i.e. the line width is dominated 
by the variance of the velocity field, which is the same for He II and H I. This is turbulent 
broadening. It has been shown recently that the velocity field v{x) shows the feature of a 
fully developed turbulence (He et al. 2006). 

To test this point, we identify the absorption lines by the similar way as that for the 
real sample. We use AUTOVP code (Dave et al. 1997) to decompose the transmitted flux 
with Gaussian proflle, and estimate the parameters of line width, column density, and the 
centroid wavelength of each line. Figure 2 shows 6(H I) vs. 6(He II) for both real data and 
simulation samples. The real data are taken from Zheng et al. (2004). Obviously, the plot 
of b(R I) vs. 6(He II) from simulated data does not follow the thermal broadening relation 
6(H I) = 26(He II), but similar to the turbulent broadening. 



5. The ratio between He II and H I Lya optical depths 
5.1. The scatter of optical depth ratio 

If the Voigt profile can be approximated by a Dirac delta function, eq.(l) yields 

7'HeIl(^) _ 1 nnciijz) 
Tmiz) 4 nmiz) 



(2) 



If atoms and ions of H I, H II, He I, He II and He III are in state of photoionization 
equilibrium, the ratio of He II to H I is (Fardal et al. 1998) 

^Hcii(^) . . ^ ^Hi 3 + ^4 / T(2:) y °^ 



nm{z) JHeii3 + Q;i VlO^'^ 
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where we assume that the UV radiation backgrounds around wavelengths wavelength c/ i/ni = 
912A and c/z/Heii = 228 A are, respectively, = Jhi('^/'^hi)~"^ and J,y = JHcii('^/^^Heii)~"'*, 
the parameters Jhi and Jneii being the specific intensities, ai and the index of their power 
laws. 

If the UV radiation background is spatially uniform, i.e. the parameters Jhi, Jneii, cti 
and a4 are constant, the ratio nudiiz) / nm{z) of eq.(3) is approximately spatially constant, 
because the temperature-dependence (T^'^^) is very weak. Thus, from eqs.(2) and (3) we 
may expect that the ratio of optical depths rHcii(-2)/THi(-2) should basically be constant, i.e. 
the scatter of ratio Tiieii{z) / Tm{z) with respect to its mean (THeii(2;)/THi(^)) has to be very 
small. 

Thus, the observed rj scatter (Kriss et al. 2001; Smcttc ct al. 2002; Zheng et al. 2004; 
Shull et al. 2004) may indicate that the ratios, Jm/JucU and (3 + a4)/{3 + ai) are not 
constant, but significantly different from pixel to pixel, or from line to line. However, eq.(2) 
is based on the assumption that the effects of thermal broadening and peculiar velocities are 
ignored. As has been shown in last section, the effects of thermal broadening and peculiar 
velocities may not be always small. Therefore, we should estimate the 77-scatter caused by 
thermal broadening and peculiar velocities. 

Since the assumption of eq. (2) is not hold, we use 77 to stand only for 4rHeii / thi {z) below. 
Figure 3 illustrates the effect of thermal broadening and a peculiar velocity field. It shows 1.) 
top panel: the 1-D distributions of the ratio 77 of optical depths given by eq.(l), 2.) second 
panel: 77 distribution from optical depths of eq.(l), but using a delta function for the Voigt 
profile; 3.) third panel: 77 given by optical depths of cq.(l) and adding, respectively, S/N=3 
and S/N=50 noises to the He II and H I transmitted flux; 4.) forth panel: temperature T 
distribution; 5.) bottom panel: mass density field of baryon matter. 

From the panel 2 of Fig. 3 one can see clearly that if the effect of thermal broadening 
is ignored, i.e. the Voigt profile is approximated by a Dirac delta function, 77 is almost 
a constant in the whole range. In this approximation, the ratio 77 does not depend on the 
fiuctuations of temperature, mass density and velocity of baryon gas, but only on the density 
ratio nHii/^Hi- The ratio nHii/nm, however, always keeps constant, even when riuii and riui 
fluctuate strongly. The rj distribution shown in panel 2 contains a signiflcant sharp decline 
at the highest temperature region. It probably results from dielectronic recombination of 
He II that is dominant at such high temperatures, and reduces the number riHen- At less 
high temperatures, the coUision ionization rate of H I is higher than He II, accordingly, there 
appear some small bumps as those visualized in panel 2. 

However, the top panel shows that, when thermal broadening is included, the fluctu- 
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ations of rj are higher at the positions with higher temperature and higher density. These 
fluctuations are due partially to the difference between the thermal velocities of He II and 
H I. It leads to the distribution of H I to be more extended than He II. The hydrodynamical 
velocity field v{x) is the same for He II and H I. It doesn't cause the change of the ratio 
nmi/nui- The effect of thermal broadening on rj scatter has also been noted by Croft et al. 
(1997), but in their samples the scatter caused by thermal broadening is not very significant. 
It is probably because more shocks are resolved on small scales in the WIGEON simulations, 
which leads to stronger temperature fluctuations (He et al. 2004). The panel 3 of Fig. 3 
shows more signiflcant scatter of rj if the noise is added. Prom Figure 1, we have seen that 
the transmitted flux of He II is substantially affected by the S/N=3 noise. 



5.2. PDFof?7 

We now examine the probability distribution function (PDF) of r], which is calculated 
pixel by pixel as ShuU et al. (2004). Figure 4 shows the PDFs for 1.) real data; 2.) simulation 
sample of the transmitted flux from eq.(l); 3.) simulation samples of the transmitted flux 
from eq.(l), and adding, respectively, S/N=3 and S/N=50 noise to the He II and H I 
transmitted flux, and 4.) simulated samples without thermal broadening, but adding noise. 
The simulation PDFs are calculated with 100 samples of 1-D transmitted flux of He II and 
H I. The error bars arc the maximum and minimum from 100 independent noise realizations. 
The simulation data arc also binned into 0.05A to match observation. 

We see that the PDF of the simulated data without adding noise has the peak at 
log?7 ~ 1.9 which is about the same as real data (Zheng et al. 2004; Shull et al. 2004). 
However, the width of PDF of the simulated data without adding noise is much less than 
the real data. The maximum scatter of rj of the simulated data is about 2{r]). The factor 2 
is just the difference between the thermal velocities of H I and He II. Though there are some 
pixels with r] > 2{r]) caused by collision ionization, it is still less than the observed scatter. 

The PDF of simulation data is significantly improved with adding noise. It looks similar 
to the observed result. The PDF of simulation data has the same peak as real data, and 
the width of PDF is also about the same as real one. Therefore, the effect of data noise is 
substantial for the scatter of 77. This effect is especially serious on the high 77 events. Samples 
without noise contamination do not have events of 77 > 300, but the samples with noise do. 
This result is about the same as Shull et al. (2004), in which events of 77 > 460 are dropped, 
because they may largely be from the uncertainty in measuring. 

Nevertheless, we see that the PDF of simulation sample with noise is still lower than 
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the real data in the range rj < 10. This result is also similar to ShuU et al. (2004). They 
found that the events of < 10 of real data show an factor 2 excess to a Monte Carlo 
calculation. Figure 4 also shows that the number of < 10 of real data is about twice as 
large as simulation sample. However, Shull et al. (2004) found that the events of < 30 
of real data arc also more than their Monte Carlo estimation by a factor 2, while Figure 3 
doesn't show so large difference. This is probably because our simulation contains the effect 
of thermal broadening and peculiar velocity field. This point can be seen from the PDF of 
r) for simulation samples without thermal broadening, but adding noise (the filled dots in 
Figure 4). It shows there are more low- 77 events without thermal broadening. In a word, our 
simulation result indicates the excess of low-?7 events in real data, but the difference between 
real and simulation data is less than the Monte Carlo estimation. 

The effect of noise can also be seen in Figure 5, which presents the relation between 77 
and thi in the range thi > 0.01. The top panel is from real data. The middle is given by 
simulation samples with adding noise of S/N=3 to He II and S/N=50 to H 1. The bottom 
shows simulation samples without noise. The rj-Tm distribution of simulation sample with 
noise is similar to the real one. The high-r/ events at low Thi area are mainly from noise. 
This is consistent with the observed result that rj is large in void area thi < 0.05 (Shull et 
al. 2004). 

Figure 6 presents the relation between 77 and column density A^(Hl) for real data and 
simulation samples with noise. It shows the correlation between r] and density of H I: 77 is 
larger for lower column density A^(Hl), and lower for higher A^(Hl). This phenomenon is 
also shown in the He II spectrum of HS1700+6416 (Reimers et al. 2005). 

6. Power spectrum 

We now compare the power spectra of H and He transmitted flux of simulation sample 
and real data. We calculate the power spectrum by the same method as Paper 1. It gives 
easy to compare the power spectra of H and He. Moreover, the real data of He is highly 
noisy, for some pixels, the S/N ratio is lower than 1, and some pixels with negative flux. 
To eliminate the effect of these pixels on the power spectrum calculation, we should use the 
algorithm of denoising or conditional-counting as Paper I (Donoho 1995; Jamkhedkar et al. 
2001). 

The power spectra of H 1 and He II Lya transmitted flux of HE2347 are shown in Figure 
7. We plot the power spectrum of H I in top panel as an indicator of the goodness of the 
simulation sample of this paper, which is produced in a box with size larger than that of 
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Paper I by a factor of 4. The parameters f — 1 and f — S mean, respectively, the threshold 
condition for denoising to be S/N > 1, and 3. Since the real data of H I Lya transmitted 
flux has high quality, the power spectrum actually is /-independent if / > 1 (Jamkhedkar 
et al. 2005). On the other hand, there are very rare pixels of the He II data with S/N > 3, 
and therefore, only a very small number of the modes is available if taking / > 3. One can 
only use / = 1 to calculate He II power spectrum of He2347. The error bars of Figure 7 is 
estimated by the maximum and minimum range of bootstrap re-sampling. 

As expected, the simulation sample of H I transmission flux is in well agreement with 
observations on all scales less than 6v = 224 km s~^. There is no discrepancy on the smallest 
scale 6v = 28 km s^^ or length scale D = 0.28 h~^ Mpc. This result is the same as Paper 
I, and therefore, the power spectrum on small scales is insensitive to the size of simulation 
box. 

The power spectrum of the He II Lya transmitted flux of HE2347, as shown in the 
bottom panel of Figure 7, is very different from simulation samples without adding noise on 
scales less than 56 km s"^. On those scales, the power spectrum of observed sample is much 
higher than the simulation results. However, the power spectrum of the noisy samples of 
He II flux gives a good fitting to the observed sample. In other words, if we could remove the 
S/N=3 noise from the real data, the power spectrum can be well fitted with the simulation 
of the LCDM model. Thus, with He II flux, we can arrive at the same conclusion as H I: 
there is no evidence for the discrepancy between observation and simulation on scales from 
1792 to 28 km s-\ or from 0.28 to 18 h"^ Mpc. 



7. Discussions and conclusions 

The H I and He II Lya transmitted flux fluctuations of QSO absorption spectrum are 
valuable to detect the fields of baryon gas and ionizing photon fleld, and to constrain models 
of the UV radiation background. With hydrodynamic simulation samples of the ACDM 
model, we made a comparison between the model-predicted statistical features and real data 
of HE2347. It includes the power spectrum, the line width of absorption features, and the 
ratio between the He II and H I optical depths. The major conclusions are as follows: 

1.) The absorption features of He II and H I basically are turbulent-broadening. It 
should be emphasized that the observed evidence for the turbulent-broadening mainly is 
given by the absorption lines in voids (Zheng et al. 2004). This is especially supported by 
the simulation samples, which shows that the turbulence behavior of the velocity field of 
IGM is on high as well as low mass density area (Kim et al. 2005). 
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2. ) The mean of the optical depth ratio (77) of real sample can be well fitted with 
simulation samples. 

3. ) The simulation samples give a very well fitting to the power spectra of both H I and 
He II transmitted flux. There is no discrepancy between the power spectra of simulated and 
observed samples on small scales. For H sample, the scales are till to 0.28 h~^ Mpc. For He, 
it is 1.1 h~^ Mpc. Therefore, there is no evidence that the power of observed sample is less 
than the standard ACDM model on small scales. 

4. ) The last but not the least is on the 77-scatter. The large 77 scatter are generally 
attributed to the inhomogeneity of the UV photon distributions. We showed, however, that 
a significant part of the T^-scatter is from the fluctuations of the temperature and velocity 
fields of IGM and data noise. There seems to be a discrepancy on low-r/ events comparing 
our simulation and the real data. The simulation sample does not contain enough low-?7 
events to fit real data. To solve this problem models of producing extra fluctuations on Lya 
transmitted flux on scales of about 1 h~^ Mpc is needed. 

The WIGEON samples are produced with a uniformly distributed UV radiation back- 
ground, and therefore, we may consider an inhomogeneous UV radiation background to be 
a possible reason of the lack of low-?] events. However, the numerical results of the inhomo- 
geneous UV radiation background shows that the fluctuations of UV radiation background 
on small scales will yield more high-77 events (200 < 77 < 350). Therefore, at least, according 
to the current calculation, the model of the UV background inhomogeneity would not be 
helpful to solve the problem of the lack of low-r^ events. 

Note that low-77 events are related to high column density iV(HI). The lack of low-rj 
events may be due to that the WIGEON samples are produced without considering the 
detailed physical processes of ionizing sources, as some mechanical and radiation feedback 
effects would not be negligible in the area of high column density, which has the high prob- 
ability containing collapsed objects. 

Any model of producing the scatter of 77 will affect the transmitted flux of both H I and 
He II. Therefore, we can use our unifled fltting of H I and He II Lya transmitted flux to set 
some constraints on models for low-77 events. 

First, if the fluctuations are Gaussian, it will play the same role as Gaussian noise, and 
increase the power of transmitted flux on scale ~ 1 Mpc. However, this paper and Paper 
I show that the power spectrum of H I Lya transmitted flux is in good consistent with 
observation on scale 1 h~^ Mpc and less. Therefore, any increase of the power of simulation 
sample will no longer be consistent with observed sample. Therefore, there is small room for 
models of adding Gaussian fluctuation. 
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Second, if the fluctuations are non-Gaussian, it should be hmited by the observed non- 
Gaussian features of the H I and He II flux. The non-Gaussian statistical features of H I 
transmitted flux, such as high order moments of the fluctuations of the flux, are sensitive 
to the addition of non-Gaussian inhomogeneity. It has been shown in Paper I that the non- 
Gaussian features of the H I transmitted flux of HE2347 can already be well fitted with the 
WIGEON samples to eight order and on scale as small as 0.28 h^^ Mpc. Therefore, there is 
also small room for models of adding non-Gaussian fluctuations. 
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Fig. 1. — Simulated samples of H I (top) and He II (second from top) transmitted flux 
without adding noise, and the same samples of H I (third from top) and He II (bottom) with 
adding noise with S/N=50 for H I and S/N=3 for He II. 
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Fig. 2. — Doppler parameter b(H I) vs. b(He II) of simulation samples (left) and real data 
of HE2347 (right). The diagonal line is expected for turbulent broadening and the dashed 
line is for b(H I) = 2b(He II). 
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Fig. 3. — 1-D distribution of rj with (top) and without (next top) thermal broadening. The 
third panel shows rj with adding S/N= 50 noise to H I and S/N=3 noise to He II flux. The 
temperature and baryon matter density are also shown. 
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Fig. 4. — PDF of Tj. 1. Real data (dot line); 2. Simulation samples with adding S/N= 
50 noise to H I and S/N=3 noise to He II transmitted flux (solid line), errorbars being the 
maximum and minimum over 100 independent noise realizations; 3. Simulation samples 
without noise (dashed line); 4. Simulation samples without thermal broadening, but adding 
noise. 
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Fig. 5. — 77 vs. optical depth thi for real data (top); simulation samples with noise (middle) 
and without noise (bottom). 
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Fig. 6. — T) vs. column density A^(HI) for simulation data (left) and real data of HE2347 
(right). 
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Fig. 7. — Power spectrum of 1.) H I transmitted flux of HE2347 with the conditional 
counting parameter / = 3 for real data (*) and simulation sample (O) (top); 2.) He II 
transmitted flux of HE2347 with the conditional counting parameter / = 1 for real data (*), 
simulation sample with adding noise of S/N=3 (A) and without noise adding (O) (bottom). 



